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1 Introduction 



The subject of three-body (and, in general, iV-body) relativistic phase space is as old as 
the hills and one might well think that all that there is to know is known about it. In 
numerical and experimental terms this is indeed true: for a long time Dalitz plots Oj 
have been routinely used in picturing data and they prove extremely helpful for picking 
out resonant intermediate states of particular spin by their preferential population of the 
plots. In the absence of any amplitude modulation by resonances or otherwise, the plots 
are at their blandest as they just represent three-body phase space. 

One of the first comprehensive references on this subject is the paper by Almgren 
In his normalization, the integral over the iV-particle phase space is defined as 



N 



I N D \p,m l ,...,m N ) = J jnd D M(P?-™. ? )0(P°)f 4Eft-P , (LI) 



where p is the total momentum. From now on, we will frequently use the notation p = y/p 1 , 
since usually it is easy to distinguish it from the cases when the four- dimensional vector p 
is meant. As a rule, we will also omit the arguments of ijy . In four dimensions (D = 4), 
we will denote In = 1$ (this is the original Almgren's notation). More details about 
integrals in other dimensions can be found in Ref. [3] and in the rest of this paper. Worth 
mentioning is that /jy is easy to work out for odd values of D, whereas considering even 
values of D brings in elliptic functions and is more difficult. 

For kinematical reasons, it is clear that the results for the integrals have no 

physical meaning if the absolute value of the momentum p is less than the sum of the 
masses. Therefore, in what follows we will imply that all results for 1^ are accompanied 
by 0{p 2 — {mi + . . . + m-Tv) 2 }, without writing this theta function explicitly. In Refs. 
integral recurrence relations for 1^ (at D — 4) were discussed. For an arbitrary dimension 
D, the generating relation can be presented as 

#>(p, mi ,...,m N ) = J ds I<S(p, A m N . R+1 , ...,m N ) J^(^ m x , . . . , m N . R ) . 

(1.2) 

Taking into account the theta functions associated with I^-r an d Ir2i, one can see that 
the actual limits of the integration variable s in Eq. (jl.2j) extend from (j2fjl R m^j to 

n-r+1 m i) ■ Another type of integral recurrence relations for I® , with respect 
to the value of D, was considered in Ref. jH]. 

The simplest example is the two-particle phase space, N = 2. In this case, the phase- 
space integral in four dimensions can be easily evaluated as 

h = j 1 \IWMM)> (L3) 

where 

\(x, y, z) = x 2 + y 2 + z 2 — 2xy — 2yz — 2zx (1.4) 
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is nothing but the well-known Kallen function [7j. 

Using Eqs. (|1.2|) and (jl.3|) (for the case D = 4, R = 1), one can obtain the following 
integral representation [3J |H] for the three-particle (N = 3) phase space: 



9 s 3 , 

7T f (IS I 



4p 



S2 



with 

Si = (mi-m 2 ) 2 , s 2 = (m 1 +m 2 ) 2 , s 3 = (p - m 3 ) 2 , s 4 = (p + m 3 ) 2 , (1.6) 

so that si < S2 < s 3 < S4. The result of the calculation of the integral (jl.5|) can be 
expressed in terms of the elliptic integrals [HI IE] (for convenience, we collect the definitions 
and relevant properties of elliptic integrals in an Appendix), 



TV 2 



^Q + {ml + m 2 2 + ml+p 2 )E{k) 



3 Vv^Tl 

+4mim2 (p — m^) 2 — (mi — m 2 ) 2 (p + m.3) 2 — m^p + m\vri2 K[k) 
+8mim2 (m 2 + m\)(p 2 + m 2 ) — 2m\m\ — 2m\p 2 IT (a 2 , A; J 

-8mim 2 (p 2 - m 2 ) 2 n (a 2 ., fc) | , (1.7) 

where we use the following notations: 

Q + = (p+mi+m 2 +m3)(p+mi-m2-m3)(p-mi+m2-m3)(p-mi-m2 + m 3 ) , 
Q- = (p-mi-m2-m3)(p-mi+m2+m 3 )(p+mi-m2+m3)(p+mi+m2-m 3 ) ,(1-8) 



/Q- 2 (P ~ m 3) 2 - ( m l + m 2) 2 2 ( m l _ m 2 



,2 



k = \ 7T 1 a i = 7 « — ? yl ' "2 = 7 ; ^"i • I 1 - 9 ) 

V Q+ [p — m 3 y - (mi - m 2 y [m\ + m 2 ) / 

We note that in Ref. jH] the notation q±± = (p±m 3 ) 2 -(mi±m2) 2 was used. In particular, 

we have Q + = q++q , Q- = <?+-<?-+, k 2 = q+~q-+/(q++q ), a 2 = q~+/q Note that 

Q± differ by the sign of p only. 

It is clear from the definition (Jl.lj) that ^3 should be a symmetrical function of the 
three masses mi, m 3 . The representation (jl.7|) in terms of elliptic integrals is however 
not explicitly symmetrical in the masses, although it must be implicitly so. One may of 
course generate a symmetrical form by averaging the unsymmetrical-looking expressions 
over the three possible permutations of m i: but this would be "cheating" since each one 
of them should be symmetrical by itself, although this is hardly transparent. 

Note that the quantities Q + and Q_ (and, therefore, the argument k) are totally sym- 
metric in mi, m.2, m.3. (In fact, they are symmetric in all four arguments p, mi, 7712,777.3.) 
Therefore, the term containing E(k) in Eq. ()1.7|) is also symmetric. The function K(k) 
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itself is also symmetric, but its coefficient is not symmetric. We also note that the product 
of Q + and Q_ produces the quantity 



^123 = Q+Q- = P 2 - {mi + m 2 + m 3 f 



V 



-mi + m 2 + m 3 ) 



p 2 — (mi — m 2 + m 3 ) 2 p 2 — (mi + m 2 — m 3 ) 



1.10) 



that occurs in recurrence relations for the sunset diagram (see, e.g., in 0[HI|). It should be 
noted that the imaginary part of the sunset diagram is proportional to the three-particle 
phase-space integral. For instance, in the notation of Ref. JHj, Im(Ti 2 3) = — 47r~ 1 I 3 . We 



also note that p® considered in Ref. |U |H| are related to 1^' as p® 



(2 



7T 



\N+D~ND A D ) 



For equal masses, mi = m 2 = m% = m, Eq. p.7|) yields 



7T 



^2 \J(p- m )(P + 3m )\ \(P 



m)(p 2 + 3m 2 )E(k eq ) - Am 2 pK(k 



1.11) 



with 



k 



oq 



(p + m) 3 (p — 3m) 
\ (p — m) 3 (p + 3m) 



1.12) 



Some other special cases of Eq. fjl . Tj) are described in Refs. |Hl HI] ■ 

This paper is devoted to a new way of exhibiting the results in an explicitly symmetrical 
manner. To do this, we will employ another integral representation for J3, in terms of 
Mandelstam variables s, t, u ^2] and the Kibble cubic <3>(s, t, u) [T3*] . In particular, we will 
show that it is convenient to present the result f)1.7|) in terms of Jacobi Z function whose 
definition is given in the Appendix. 



2 Phase space integrals 

As an illustration, let us demonstrate how the connection with the Dalitz picture can be 
derived directly from the definition (jl.lj) . The D-dimensional vector p can be presented 
as (p°,p), where p is the (D — l)-dimensional Euclidean vector of space components. 
Without loss of generality, we can work in the center-of-mass frame, p = (p°,0). Using 
the integral representation 

S (YlPi-p) = (2^d / dD °° expji^(pix) -i(px)| , (2.1) 
with (px) = p°x°, we get 

4? } = J^y5 J d D x e~^° I ft / d D Pl 6(p 2 - m 2 ) 6(p°) } . (2.2) 
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(Similar method was used in Ref. [H].) Integrating over (D — 1) -dimensional angles of p t 
we get 



(2.3) 

with pi = |pj| and £ = |x|. In the four- dimensional case the Bessel function reduces 
to an elementary function, Ji/ 2 (pi^) = [2/(7rpi£)p 2 sin(pj£). We note an analogy with 
the calculation of Feynman integrals in the coordinate space |15| I22j. when each massive 
propagator yields a (modified) Bessel function. 

Let us consider, for example, the two-particle phase space. Then, the integration over 
£ gives us 

oo 

J £d£ J v ( Pl £) J,(p 2 = 25{p{- pi) , 
o 

with z/ = (£) — 3)/2, so that we can put p\ = p 2 = p, whereas the integration over x° 



yields another delta function, 5 yp — yp 2 + mf — y p 2 + wilj in the center-of-mass frame. 
The resulting integral 

(D-l)/2 °? n D - 2 Hn / \ 

can be easily evaluated, yielding (see, e.g., in Ref. 0) 



j(D) 



(D-1J/2 



(2p)^r(^i 



(D-3J/2 



A(p 2 ,m 2 ,m 2 ) , (2.5) 



where A is defined in Eq. fjl.4|) . For D = 4, Eq. ()2.5|) reduces to the well-known an- 
swer fjl.3|) . 

For the three-particle phase-space integral we get 

D 2 V 7T /■ /• ip O x o ] r /■ Pj dpi . o v ^f^? 

1 V 2 J S -oo 1=1 + * 

(2.6) 

Here we can integrate over £, using (see Ref. [TH]) 



T^=i JvKPxK) Jv\p2K) MpsS) , a 

4 xv 2 r[v + ±)(8pip 2 p 3 y 



(with v = (D — 3)/2), where A is the Kallen function (jl.4|) . In fact, in our case, when all 
P* > 0, 

9{-\{pl pi pi)} = 6( Pl + p 2 - p 3 ) 9(p 2 + p 3 - pi) 9(p 3 + Pi - p 2 ) , (2.8) 
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i.e., it equals 1 when one can compose a triangle with sides pi, p 2 , P3, and gives otherwise 
(cf. Eq. (11) of Ref. [T7j). 



Introducing notation <jj = J p 2 + mf and integrating over x° (getting a 5 function) we 
arrive at 



7T 



D-2 



oo oo oo 



4 D) = r p _ 2 ^ / / / 'W da :j M p - ct i - (T-, - o- :s ) 



X 



mi rrt2 m,3 

(Z?-4)/2 



\/ 2 2 2 2 2 2\1 ^ a ( \/ 2 2 2 2 2 2\1 /o n\ 

A((T 1 -m 1 ,a 2 -m 2 ,cr3-m 3 )J {-\{cr 1 -m 1 , a 2 -m 2 , a 3 -m 3 ) j .(2.9) 



In four dimensions the factor [— A]^~ 4 ^ 2 disappears and, geometrically, we need to cal- 
culate a closed area on the plane <i\ + o 2 + a 3 = p° = p, with the boundary of the figure 
described by 

A((7i - m\, o\ - m 2 , a 2 - ml) = 0, <7i + a 2 + cr 3 = p . ( 2 -10) 
Furthermore, introducing Mandelstam-type variables 

Q Q OO OO /\ 

s=p + m 3 — 2p<T 3 , t = p +m l — 2p<Ji, u = p + m 2 — 2pa 2 (2.11) 
satisfying 

s + t + u = m\ + m\ + ml + p 2 = w , (2-12) 

one arrives at another integral representation (the limits of integration are discussed 
below) , 



J J J dsdtdu5(s + t + u-w ) [$(s,t,M)] (D " 4)/2 0{$(s,t,M)} , 

(2.13) 

— A {\(s,m 2 3 ,p 2 ), A(t,m 2 ,p 2 ), X(u,m 2 2 ,p 2 

)} (2-14) 



where 

•(«.«.») Vif - 

can also be written in a more familiar Kibble cubic form |13j . 

$(s, t, u) = stu — s(m 2 m 2 + p 2 m 2 ) — i(m 2 m 2 + p 2 m\) — u(mlm\ + p 2 m 2 ) 

,o/ 222i222i222i22 2\ /oicr\ 

+2(m l m 2 m 3 +p m l m 2 +p m 2 m 3 +p m 3 m 1 ) , (2.15) 
provided that the condition ()2.12|) is satisfied. In particular, in four dimensions we have 

I =^ 
3 4p 2 



ds dt dw <J(s + t + u - w ) 9 {$(s, t,u)} . (2.16) 



According to the definition 1)2.11)1 in terms of cr,, one can see that the maximal values 
of s, t and u (corresponding to the upper limits of integration in Eqs. (j2.13j) and (|2.16jl ) 

are s max = (p — m 3 ) 2 , t max = (p — mi) 2 and M max = (p — m 2 ) 2 . To define the minimal 
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values of s, t and u, the familiar Dalitz-Kibble picture given in Fig. 1 is useful. Due to 
the condition ()2.12|) . the region of integration is restricted by a triangle, with s > s m ; n = 
(mi + m 2) 2 , t > t m in = (iti 2 + m 3 ) 2 and u > u min = (m 1 + m 3 ) 2 . Moreover, due to the 
theta function 8 {Q(s,t,u)} the region of integration is in fact restricted by the interior 
of the cubic curve Q(s,t,u) = 0, see Fig. 1. Within that area, the Mandelstam variables 
s, t and u take their minimal values in points O s , Ot and O u , respectively, whereas their 
maximal values correspond to the points P s , Pt and P u . (The dashed triangles will be 
discussed below in Section 4.) 

The function Q(s,t,u) has a maximum within the region of integration. For equal 
masses, the maximal value $ max = ^ p 2 (p 2 — 9m 2 ) 2 occurs at s = t = u = | (p 2 + 3m 2 ). 
For the general unequal masses, one needs to solve a fourth-order algebraic equation to 
find the position of the maximum. 

We note that the representation (j2.14|) can be extracted from Eq. (5.39) of [T%] . using 
symmetry properties. Our $(s,t,w) corresponds to — G(s, t,p 2 , m 2 ,, m 2 , m 2 ), in the nota- 
tions of [IB]. The G-function is symmetric with respect to the permutations of three pairs 
of arguments, (s,t), (p 2 ,m 2 ) and (m 2 ,m 2 ). Although the authors presume from their 
Eq. (5.39) that "from a practical point of view this identity is not very useful", we found 
that its symmetric form is certainly helpful in understanding the structural properties of 
phase space integrals. 
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3 Geometrical interpretation 



Let us introduce 



t 



mr, 



m% 



u 



ml 



m% 



Cl2 



C23 



Cl3 



2mim2 ' 2m 2 m% 2mim 3 

Then, the function $(s, t,u) can be presented as a Gram determinant, 



(3.1) 



whereas the 5 function becomes 
5(s + t + u 



A 2 2 2 

4m 1 m 2 m 3 



1 C12 C13 
Cl2 1 C 23 
C13 C 23 1 



(3.2) 



m 



2 -m 3 -p 

2 1 ™2 1 2 



5 mj + m 2 + m 3 + 2mim 2 ci2 + 2m 2 m 3 C2 3 + 2mim 3 Ci 3 — p ) . (3.3) 



In this way, we get 



2VI" o r) o 

— —m,m 9 mn 
p 



-12 





/ 


1 C12 C13 




dci 3 dc 23 




Cl2 1 C 23 


) 




V 


C13 C23 1 





x5 (ml + ml + ml + 2m l m 2 Ci 2 + 2m 2 m 3 c 2 3 + 2m 1 m 3 Ci 3 - p 2 ^ 



(3.4) 



where the integration extends over Cji > 1. 

If one were to interpret Cji as the cosines of the angles between the mj and mi sides of a 
vertex of a parallelepiped (formed by mi, m 2 and m 3 , see Fig. 2), then all these quantities 
would have a straightforward geometrical interpretation. Namely, Q(s,t,u) would be 
4{volume of parallelepiped} 2 , whereas the 5 function would tell us that the "principal" 
diagonal of this parallelepiped should be equal to p. In this case, the quantities y^i, y/t 
and y/u could be identified as the diagonals of the faces of the parallelepiped, see in Fig. 2. 
Moreover, 

p 2 + m? — t p 2 + m\ — u . p 2 + m 2 , — s . . 

, and 2 (3.5) 

2pm\ 2pm 2 2pm 3 

could be understood as cosines of the angles between the diagonal p and the mj sides of 
the parallelepiped. In other words, these are the angles between p and rrii in triangles 
with sides (p, mi, v^), (p,m2, y/u) and (p, m 3 , a/s), respectively. 

Using this geometrical picture, we can mention a rather interesting geometrical mean- 
ing of Eq. (j2.14J) . Namely, it tells us that the volume of the parallelepiped is (8/p) times 
the area of triangle whose sides are given by the areas of triangles formed out of the 
principal diagonal p, one of the face diagonals (\/s, \/t or y/u), and the appropriate m 3 , 
mi or m 2 side. 

However, when we are above the threshold, p 2 > (mi + m 2 + m 3 ) 2 , the quantities Cji 
exceed one and therefore the expressions should be understood in the sense of analytic 



8 



Figure 2: The parallelepiped interpretation 



continuation, i.e. as hyperbolic cosines. The same is valid for the triangles (p, m 3 , y/s), etc.: 
they should also be understood in the sense of analytic continuation, since p > m 3 + y/s, 
etc. Therefore, the quantities Oijvrii should also be understood as hyperbolic cosines, 

whereas \j of — mf / are hyperbolic sines. 

Nevertheless, in the region below the threshold (which we need, for instance, to de- 
scribe the real part of the sunset diagram), this geometrical picture can have direct mean- 
ing, generalizing the picture we had for the one- loop two-point function |19j . 

4 Kibble cubic characteristics 

Suppose 

(s ,t ,w - s - t ), (so, w - s - Mo, Mo), (w - t - u ,t , u ) (4.1) 
all are the roots of the equation $>(s,t,u) = 0. Then, we can present $>(s,t,u) as 

$(s, t, u) = stu — stoUo — Sotuo — sotou + 2sotoUo ■ (4.2) 
Furthermore, if we shift the Mandelstam variables as 

s = so + s' , t = to + t' , u = uo + u' , (4.3) 
subject to the condition 

s' + t' + u' = wo - s - t - uo = w' , (4.4) 
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then 



$(s, t, u) s't'u' + s t'u' + s't u' + s't'uo . 
Using equation Eq. ()4.2|) and denning 



Ctu = 



/to «o 



C s t 



's t 



C.C7/, 



S U 



tu V st V su 

we arrive at another Gram determinant representation for $(s, t,u) (cf. Eq. ()3.2|) ' 



(4.5) 



(4.6) 



$(s, t, u) = stu 



1 Qu c s t 

Ctu 1 C su 
Cst C sn 1 



(4.7) 



There are (at least) two sets of solutions (|4.1|) that can be described as 

A 2 A 3 



A X A 2 

S — ' A j 



tn 



A X A 3 
u o = — — , 



M A x ' A 2 

so that Eq. (l4~2l yields 

$(s, t, w) = stu - A\t - A\u - A\s + 2A X A 2 A 3 . 

The first set of solutions corresponds to 

A\ = prrii + m 2 m 3 , A 2 = pm 2 + m 3 m\ , A 3 = pm 3 + m\m 2 

For this set, we have 

m\m 2 m 3 p Q- 



w = w - s - t - Uo 



Ctu 



pm 3 + mi 777.2 



A±A 2 A 3 

pm 2 + 777.1777.3 



]?777i + 777 2 7773 



-st 



(4.8) 

(4.9) 
(4.10) 

(4.11) 
(4.12) 



Note that if we change p — > — p in Eq. f)4. lUjl . this would also be a solution, which would 
correspond to a "non-physical" branch of the Kibble cubic. 
The second set of solutions corresponds to 

A\ = \{p 2 + m\ — m\ — 777,3) j A 2 = \(p 2 — m\ + m\ — m\) , A 3 = \(p 2 — m\—m\ + m\) . 

' (4-13) 

(4.14) 
(4.15) 



For this set, we get 

w = w - s - t - m 



D 



123 



16AiA 2 A 3 I6A1A2A3' 



p—mf— 777,2+7773 p^ — m\ + 7773— 7773 



2Vtu 



2^sl 



p +1711—1712 — m 3 
2\/su 
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It should be noted that the value of s corresponding to the second set satisfies 



X(s ,ml,m 2 2 ) = \(s ,p 2 ,ml) , (4.16) 

i.e., the areas (or their analytical continuations) of triangles with sides (^ l /so, m\, m 2 ) and 
(y/so, p, m.3) are equal. Analogously, 

X(t ,ml,ml) = X(tQ,p 2 ,m\) , A^o.m^mj) = X(u ,p 2 ,ml) . (4.17) 

Moreover, one can get the direct geometrical interpretation of the quantities ()4.15j) through 
the familiar parallelepiped picture given in Fig. 2. Namely, Ct u is nothing but the cosine 
between the face diagonals \Ji and \fu. Accordingly, c su is the cosine of the angle between 
y/s and \fu diagonals, whilst c st is the cosine of the angle between y/s and \ft diagonals. 
If we construct a tetrahedron using the y/s, \/t and y/u diagonals then, according to 
Eq. ()4.7j) . <3>(s, t,u) would represent 36 times its volume squared. 

In the Dalitz-Kibble plot shown in Fig. 1 we connect the points (J4.1J) for each of the 
two sets by dashed lines, introducing subscripts [1] and [2] for the first and the second set, 
respectively. The two resulting "dashed" triangles indicate that the two sets of solutions 
are complementary to each other. Namely, the boundary of the Dalitz plot confines the 
products tu, st and su as follows: 

(pm 3 + mim 2 ) 2 <tu< \{p 2 — m\ — m 2 , + ml) 2 , 



(777772 + mim^) 2 < st < j(p 2 — ml + m 



2 ™2\2 



2 



m, 



(pmi + m2fns) 2 < su < \{p 2 + m\ — m\ — m\) 2 , (4-18) 

or, equivalently, 

(*0«o)[l] < tU < {t U )[2} , (s *o)[l] < St < (S t )[2] , (s Mo)[l] < SU < (s U )[2] • (4-19) 

In other words, the first and the second sets yield, respectively, the minimal and the 
maximal values of tu, st and su. 

Let us consider the corresponding values of the "cosines" c su , c st and c tu . For the first 
set, c su , c st and c tu would vary between 1 and cosipi (i = 1,2,3), respectively, where 

2(pmi — 17121713) 2{pm,2 — 777.3777.1) 2Qom 3 — 771^2) 

COS^i = ^- 2 2 2> c OS<^2 = ^ 2 2' COS (f 3 = — r f . 

p iJ rV(l\— 1712— m 3 V mf + m^ — 7773 77 z — 777^ — 7772 + 7773 

(4.20) 

For the second set, c su , c st and c tu would vary between 1 and 1/ cosy?,. This means that 
we need to understand them in the sense of analytic continuation. 

The angles <^ will be very important below. Their sines can be presented as 

VQ+ ■ VQ+ ■ VQ+ 

sm^i = — — = J, sm^ 2 =^ s- — o 2> S1I1( ^3 



p 2 +mi — 777-2 — m 3 J9 2 — 777i + 7772— 7773 p 2 — m( — 1712 + 171^ 

' (4-21) 
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It is interesting that the corresponding Gram determinant can be factorized as 



1 — COS <y? 3 — cos (f 2 
cos (p 3 1 — cos Ifi 
cos <p 2 — COS (fi 1 



— sin 2 y3i sin 2 yj 2 sin 2 ^3 



(4.22) 



Eq. ()4.22|) can be used to express k in terms of </?j. We also note that 



tan 



(p — m 3 ) 2 — {nil — 1^2)' 
\ (p + m 3 ) 2 - (mi + m 2 ) 



2 ' 



(4.23) 



and similarly for ipi and y? 2 - In particular, one can see that at the threshold, p = mi + 
m 2 + m 3 , the angles <fi are related to the angles 6^ from Eq. (20) of Ref. [20] (see also in 
Ref. [10 J as(pi = n-29i, and 



We can also consider associated angles ipi, such that 



7T . 



sin i/ji = k sin <y2j, cos ^ = \ll — k 2 sin 2 . 



Explicitly, we get 

sin ip 3 



9 2 2 1 2 ' 

£r — — m 2 + m 3 



cos ^ 3 



2(pm 3 + mim 2 ) 

9 2 2 I 2 ' 
p z — rrii ~ m 2 + m 3 



etc. For these angles, we get 



1 cos ^3 cos -02 
cos ips 1 cos -01 

COS -02 COS -01 1 



— sin 2 ipi sin 2 -0 2 sin 2 ip 3 



(4.24) 



(4.25) 



(4.26) 



(4.27) 



5 A naturally symmetric representation 

Using the representation ()4.2|) for $(s,i, u), in terms of so, t Q and uq, the three-body 
phase-space integral can be written as 

I3 = ds dt du 5(s + t + u — wo) 8(stu — st uo — s tu — sot u + 2sot uo) , (5.1) 

with wq = p 2 + m 2 + m 2 + m 2 . Integrating over u yields 



7l~ 

Ap 2 



ds dt 9{(st — sot )(wo — s — t) — st uo — Sotuo + 2sotoWo} • (5.2) 
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Then, integrating over t, we basically obtain the difference between the roots of the 
quadratic argument of the 9 function, which is 



- \J s 4 - 2w s 3 + (w% + 2s t + 2s u - 4t M ) s 2 - 2 (w t + w u - 4:U t ) s s + si (t -u ) 2 . 

It is easy to check that for both sets of (so,to,Uo) the square root takes the familiar 
form (jl.5|) . which yields the non-symmetric result f)1.7|) in terms of elliptic integrals. 

Starting from the representation (|2.13|) . one can easily generalize the result ()1.5|) to 
the Z)-dimensional 

^ = W »S(-») 7?^ [(-*)(-*)(•.- ■')(«- -r- 3 " 2 ■ (5-3) 

with Si given in Eq. (|1.6p . Another way to derive the representation ()5.3|) is to use the 
recurrence relation (|1.2jh 

4 D) = y ds if ) (p, vG,m 3 ) 4 D) (v^,^i,^ 2 ) , (5.4) 

and substitite the result ()2.5j) for 1% . The result ()5.3)1 corresponds to Eq. (9) of Ref. jl]. 
[We note that the overall factor on the r.h.s. of Eq. (9) of Ref. [I] should be corrected: 
(32tt) 2 - 2£ should be changed into ±(167r) 2 ~^, with t = D/2]. 

Using representation ([5.3)1 . it is easy to see (just substituting s = x 2 ) that all odd- 
dimensional phase-space integrals can be expressed in terms of polynomial functions (see, 
e.g., in Refs. [HI 121 123), 

4 3) = Trip ~ m i ~ m 2 - m 3 ) , (5.5) 
2p 



J ' 5) = ^f^~ mi ~ m2 ' m ^ 



~(p — mi — m 2 — m 3 ) 4 + (mi + m 2 + m 3 )p 3 



—2{m 1 + m 2 + m 3 )p + ym l + m 2 + m 3 )p + I2mim 2 m 3 p 

— (m! + m 2 +m 3 ) (m! +m 2 ) (m 2 +m 3 ) (m 3 +m!) + 4mim2m3(m 1 +m 2 +m 3 ) 



, (5.6) 



etc., which are explicitly symmetric in the masses mj. However, the results in even 
dimensions appear to be less trivial. 

It is instructive to consider the two-dimensional case, D = 2. Then, the integral ()5.3|) 
yields just the elliptic integral K(k), 



r(2) _ 

S2 



ds 



Is = I , = -jfr K(k) . (5.7) 
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This is of course explicitly symmetric in the masses without further ado. On the other 
hand, using the 5 function in Eq. (|5.1|) . we can insert 1 = (s + t + u)/wq in the integrand, 
and then consider the three resulting terms (with s, t and u) separately. In this way, we 
arrive at an alternative expression, 

4 2) = -1 7 . sds + J . tdt 



W °U ^J(s-s 1 )(s-s 2 )(s 3 -s)(s 4 -s) + l J(t- tl )(t-t 2 )(t 3 -t)(U-t) 



V-3 . 

u du 



+ / / > (5-* 

J {u — Ux){u—U-2)(Uz — u){u±—u) ) 



(1 2 



where the roots U and Ui can be obtained from s« given in Eq. (jl.fij) by proper permutation 
of the masses mj. Each of the integrals involved in Eq. (|5.8j) can be expressed in terms of 
Jacobian Z function (see in Appendix). For example, 



J Sds = . Siny . 3 K(k) - K(k) Z(<p 3 , k) , (5.9) 



where (fi are nothing but the three angles defined in Eqs. (|4.2()jl - (|4.21j) . Comparing the 
resulting expression with the original result ()5.7|) we obtain a very useful relation between 
the three Z(<pi, k) functions, 

Z(<pi, k) + Z((f2, k) + Z({ps, k) = k 2 sin^i siny9 2 sin<^ 3 . (5.10) 

Let us now consider the four-dimensional integral 1% = 1%, namely, its representa- 
tion p.5|) . A useful observation is that the result would be simpler if we managed to 
get rid of s in the denominator. In particular, it would contain just one elliptic integral 
II, rather than two. How are we to eliminate s in the denominator? Again, using the 5 
function in Eq. (|5.1J1 . we can insert 1 = (s + t + u)/wq in the integrand, and then consider 
the three resulting terms (with s, t and u) separately. For the term with s, we perform the 
t and u integrations and arrive at the same integral as in Eq. (|1.5|) . but without s in the 
denominator. In two other integrals, we just integrate in a different order, leaving as the 
last one the t or u integration, respectively. In this way, we obtain for the integral (|5.1|) 



TV 2 



Ap 2 w 



ds dt du (s+t+u) 5(s+t+u — w ) 9(stu — st u — s tu — Sot u + 2s t u ) 



TV 2 



S3 *3 



Ap 2 wo 



ds ^/(s-s 1 )(s-s 2 ){s 3 -s){s A -s) + dt ^(t-hXt-hXh-t^U-t) 



■S2 *2 



+ / du \J(u - Mi) (u - u 2 ) (u 3 - u) (u 4 - u) \ , (5.11) 



U2 
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where, as before, the roots £j and u,i can be obtained from Sj given in Eq. by 
permutation of the masses. 

Using the formulae given in Ref. [21] along with Eqs. (jA.8|) and (jA.ll|) . the s-integral 
in Eq. ()5.11|) can be calculated in terms of a Jacobian Z function (see Appendix), 



A' 2 



!Q + { 2{p 2 m\ - mlml)K(k)^^- + ±Q.K(k) 



+ 



6 



(p 2 -m 2 1 ~m 2 + ml) 2 + 8(p 2 m 2 3 + m\m 2 2 ) [E(k) - K(k)} \ , (5.12) 



where f% is one of the three angles defined in Eqs. ()4.20|) - ()4.21|) . 

Collecting the results for all three integrals and using the relation ()5.10|) . we arrive at 
the symmetric result 



7T 



^ 3 = o~2 i V<3+(P +mi + m 2 2 + m 2 3 ) [E(k) - K(k)} 



+Q+K{k) 



Z(<Pi,k) Z((p 2 ,k) Z(f 3 ,k) 

. 9 ~r . 9 ~r . 9 



sin <f i 



sin if 2 



sin if 3 



(5.13) 



This symmetric result can also be presented in terms of the elliptic integrals II, using (see 
in [21) 

K{k) Z(ipi, k) = cot <fi \Jl^- k 2 sin 2 ifi [u(k 2 sin 2 <f h k) - K(k)] . (5.14) 

In principle, one can also derive the result ()5.13|) directly from the non-symmetric 
representation (|1.7jl (see in Ref. [22]), in a tedious way relying on the use of several 
relations collected in the Appendix, including the addition formula ()A.9j) for Jacobi Z 
functions. 

Worth noting is that in a similar way one can obtain results for higher even dimensions 
D. For instance, in six dimensions we get 



r(6) 



7T 



Q 



1/2 



lUp 4 1 20 



[E(k)-K(k)\ 192(p s +ml + ml + ml) - ll2(p 4 + m{ + m 2 +m 



4\2 
3) 



— 6(p 2 + m 2 +m 2 + m 3 ) 4: — 156(p 6 +ml+m 2 +ml) (p 2 +m 2 +m 2 +m 3 i ) 
+83(p 4 +m 4 +m 2 +^3) {p 2 +m\ +m 2 +m 2 ) 2 



1 

40 



+ 1 ^Q-Q 1 / 2 K{k) \3{p 2 + mj + m 2 2 + m 2 ) 2 - 16(p 4 + m 4 + m 4 + m 4 ) 



+ 3 QfK(k) 



4 sin<£>i sin if 2 sin<£>3 



Z{fi,k) Z(if 2 ,k) Z(f 3 ,k) 

. 9 ~r . 9 ~r . 9 



sm if 1 sin if 2 
15 



sm f 3 



-qQ+(P +mi + m 2 2 + mi)K{k) 



Z(ipx, k) Z(<p 2 ,k) Z(<p 3 ,k) 

+ . A ~T . A 



sin 4 ifx 



sin ip 2 



sin y? 3 



. (5.15) 



As an alternative way to obtain results for higher values of D, the approach of the paper [9J 
may be used. 

In the equal-mass case, 



<Pl = ¥2 = V?3 = <Pe q , Sin V?eq 

Here, using Eq. (|5.10|) we get 



(jy — m)(p + 3m) 



p + m 



COS y?eq 



2m 
p + m 



Z((f eq ,k eq ) = - k 2 sinVeq 



(5.16) 



(5.17) 



with k eq defined in Eq. (11.12)1 . In this way, we reproduce Eq. (jl.ll)) . whereas for D = 6 
Eq. (EHo) yields 



7T 



. ? v / (p-m) 3 (p + 3m){ (p 4 - 9m 4 )(p 4 - 42p 2 m 2 + 9m 4 ) - K(k cq )} 

+{p + m) 3 (p - 3m)(p 4 - 36p 2 m 2 + 27m A )K(k cq )} . (5.18) 



We note that Eq. (J5.17)) yields a reduction formula of Z(ip, k), for a special case when 

k 



\/l — 2 cosy? 



sin<yj(l — cosy?) 

Another interesting limit corresponds to the case when one of the masses vanishes (for 
example, m 3 — > 0). This corresponds to the case k — > 1, when -E(/c) is finite (£"(1) = 1) 
whereas K(k) develops logarithmic singularity. At m 3 = 0, cosy? 3 < and y? 3 > tt/2, so 
that we need to use Eq. (|A.7|) . Using equations listed in Ref. we get 



lim {K(k) [±Z{<p, k) - sin p]} = -\ In (l±^ 
k-*i 2 \ 1 — sin (p 



(5.19) 



where plus or minus should be used for p < tc/2 or p > tc/2, respectively. Let us consider 
Eq. (j5.13j) . Using Eqs. ()5.19j) and (|4.21|) we see that singular terms containing K{k) 
cancel, and we arrive at the following result: 



lim J 3 

TO3— »0 



7T 



57r 



l Q + {p 2 + m 2 + m\) H — (p 2 — m 2 — m 2 ) 2 In 



2^ - /;; T - >>>2> 

■3 1 ™2 ^,2 



'p 2 — m 2 — 777-2 + \[Q~ 



9 2 2 

7T — mf — TOj 



2 , ^2 ^2 ln 



--(p +m x -m 2 



•-(p 2 — m 2 + m 2 ,) 2 ln 



p 2 + mf - 777-2 + vQ^ 

y 77 2 + 777 2 — m 2 — 

2 - m\ + m 2 2 + y/Ql 



P 



y p 2 — m 2 + m 2 
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(5.20) 



where Q + = X(p 2 , m 2 , m 2 ,) in this limit. It is easy to check that this expression is equivalent 
to known results (see, e.g., in Refs. j3J|H]). The advantage of our approach is that the 
symmetry with respect to any of the remaining masses is always explicit, whereas non- 
symmetric expressions like Eq. (jl.7j) lead to the answers which are not explicitly symmetric 
(cf. Eq. (57) of Ref. 0). 

6 Conclusion 

We have considered several representations for the three-particle phase space, exploring 
their symmetry properties and geometrical meaning. It was shown that the angles (pi 
defined in Eqs. (j4.2Uj) - (j4.21j) are convenient to describe the results for the three-particle 
phase-space integral 1%. In terms of the Jacobian Z function (related to the elliptic 
integral II through Eq. (j5.14j0 . the result for J 3 in four dimensions is given in Eq. (j5.13|) . 
It is very compact and explicitly symmetric with respect to all masses m«. Note that 
the three zeta functions Z(ipi, k) are connected through the relation (j5.10j) . This relation 
can be obtained by comparing the representation ()5.7|) for two-dimensional integral 1^ 
with another representation obtained by using the delta function properties. We have 
also considered the six-dimensional case. The result for I3 is given in Eq. (|5.15jl . also 
expressed in terms of Z(ipi, k). 
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Appendix: Elliptic integrals 

The normal elliptic integrals of the first and second kind are defined as 






(A.2) 



At ip = 7r/2 we get the complete elliptic integrals, 



1 



dt 



( 



1 1 



) 



K(k) 




I 



(l-t 2 )(l-Pt 2 ) 





(A.3) 












(A.4) 
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n( c , k) 



At 



(l-ct 2 )J(l-t 2 )(l-k 2 t 2 ) 



7T 

2 



F, (§; 1, |; 1| c, fc : 



(A.5) 



where Fi is the Appell hypergeometric function of two arguments. 
The Jacobian zeta function, Z(f3, k), is defined through 

K(k) Z(J3, k) = K(k) E(fi, k) - E(k) F([3, k). 



(A.6) 



We will assume that < k < 1. (In the limit k — *■ 1, F(/3,k) and are singular.) 

From the definition (|A.6|) it is obvious that Z (^,k) = 0. Moreover, using symmetry 
properties of E({3, k) and F({3, k), we get 



Z { f + 5, fc 



(A.7) 



To represent the elliptic functions H(a{, k) occurring in Eq. f)l .7|) in terms of Z func- 
tions, we can use 



Tl(a 2 ,k) =K{k) + 



a, K(k) Z(/3i,k) 



(A.8) 



^{l-a}){k 2 -aj) 

with Pi = arcsin(«i/A;). Eq. (|A.8|) corresponds to case III on p. 229 of [21], when < 
a 2 < k 2 . 

The following addition formulae from Ref. [21] (p. 34, Eq. (142.01)) are needed: 



Z(Pi, k) ± Z((3 2 , k) = Z((p±, k) ± k 2 sin (3\ sin /3 2 sin (p±, 



(A.9) 



where the angles 



2 arctan 



sin /3i \Jl — k 2 sin 2 /3 2 ± sin /3 2 \j\ — k 2 sin 2 ft 

cos ft + cos ft 



(A.10) 



In fact, the angles y?_ and <£>+ correspond to the angles (fx and <£> 2 (see Eqs. (|4.2()jl - (|4.21|l ). 
respectively. Moreover, using the same addition formula (|A.9J) with ft j2 substituted by 
<fi,2, we get the symmetric connection ()5.10|) between Z((fi,k) (i = 1,2,3). 

To derive the result given in Eq. ([5. 12)1 . one can use the integral tables of Ref. 
along with the following relation: 



2Z{/3 1 ,k) = -Z(<p 3 ,k) 



2k 2 sin 3 ft cos fty 1 - k 2 sin 2 ft 
l-k 2 sin 4 ft 



(A.ll) 



It corresponds to the last two lines of Eq. (141.01) on p. 33 of Ref. |24j . where ip <->• 77 — ^3. 
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